Quantum Films Adsorbed on Graphite: Third and Fourth Hehum Layers 
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Using a path-integral Monte Carlo method for simulating superfluid quantum films, we investigate 
helium layers adsorbed on a substrate consisting of graphite plus two solid helium layers. The solid 
helium layers are modeled first as inert, with paths frozen at equilibrated positions, and then as 
active, with second layer atoms included in the Monte Carlo updating. In both cases, we observe the 
formation of as many as three well defined additional layers above the first two, and determine the 
layer promotion density by calculating the density profile and through a calculation of the chemical 
potential. For liquid layers adsorbed onto the inert solids, we find self-bound liquid phases in both 
the third and fourth layers and determine the equilibrium density. In the third layer at coverages 
below equilibrium, we find liquid droplets and a metastable uniform liquid phase, and determine 
the spinodal point that separates these regions. The above phases and their coverage ranges are in 
good agreement with several experiments. The superfiuid density as a function of coverage is also 
calculated, and it is observed to change only weakly around the promotion density. For coverages 
above the beginning of fourth layer promotion, we observe continued increase in the third layer 
density. We note that the third layer density increase is perhaps enough to cause solidification 
in this layer, which would explain heat capacity peaks observed experimentally for fourth layer 
coverages and would provide a simple explanation for the plateaus seen in the superfluid coverage. 
For helium adsorbed on an active second layer, we observe that a self-bound liquid phase occurs 
in the third layer, and we determine the equilibrium density and spinodal point, which remain in 
agreement with experiment. We find that promotion to both the third and fourth layers is signaled 
by a change in the density dependence of the chemical potential. We further observe the increase in 
the second layer density with increasing total coverage. The coverage dependence of the superfluid 
density is calculated, and a pronounced drop is seen at high third layer coverages, as has also been 
observed experimentally. 
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I. INTRODUCTION 

Helium films adsorbed on graphite exhibit a number 
of phases and have proven to be a rich source for both 
experimental and theoretical studies of two-dimensional 
phenomena. The graphite substrate is ordered on atomic 
length scales and offers a potential well for helium that 
is relatively strong for physical adsorption but is short 
ranged perpendicular to the substrate. As a result, a 
number of distinct, atomically thin layers occur, each 
with its own phase diagram. Near the graphite surface, 
the layers tend to solidify, with both commensurate and 
incommensurate solids occuring in the first two layers. In 
both layers, this solidification occurs before promotion to 
the next layer. The second layec exhibits a coverage re- 
gion with superfluidity as well,Eru while the first layer 
apparently favors the formatipa of solid clusters over liq- 
uid droplets ^trjow densities jETlI a Ithough there is debate 
on this issucaiJ A general review of physically adsorbed 
films such as helium on graphite nari be found in the 
book of Bruch, Cole, and Zaremba.E3 In this paper, we 
will focus on the liquid third and fourth layers. 

Detailed information on the film structure at low tem- 
peratures for the third and higher layers has come-fmm a 
number of experimpats, including heat capacity,Eja tcp^ 
sional oscillator ,Wll3 and third sound measurements. 



Ranges in which the heat capacity depends linearly on 
the coverage suggest that gas-liquid coexistence regions 
exist in the third and fourth layers. The transition from 
liquid droplets to a uniform liquid phase in each layer jia 
signaled by a peak in the isothermal compressibility.ES 
Unlike the first and second layers, these higher layers do 
not solidify before layer promotion, since torsional oscil- 
lator measurements detect superfluidity for all coverages 
beginning at intermediate third layer densities. This ap- 
parently ruleS|Xi.ut earlier suggestions that the third layer 
may solidify,Llllj although it still may be possible for solid- 
ificatino to take place under compression of higher layers. 

Perhaps the most unusual feature of the higher layers 
is the step-like behavior of superfluidity with increasing 
density.EJI3 The superfluid coverage in layered films will 
not grow continuously because of the layering transitions. 
Plateaus in the superfluid density immediately after layer 
promotion are expected, and have been observed, when 
the particles are in the droplet region. These plateaus 
occur because the droplets in the new layer lack the cop^ 
nectivity to exhibit superflow across the entire surface.llil 
The interesting observation is that the plateaus actually 
begin before layer promotion for the third through sixth 
layers. At 500 mK, superfluidity even exhibits a decrease 
with increasing coverage near the promotion to the fourth 
layer. These effects have been discussed in the context 
of the Bose-Hubbard model.tj The suggestion is that the 
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plateaus are produced by the increased localization of 
particles in the dense liquid. 

Theoretical tools applied in the study of quantum 
films on realistically treated substrates include the-hii^ 
pernetted chain Eulcr-LagrangaiHNC-EL) mcthod,EJo 
densits2-| ftuiCjtional theoryJ£3~E3 and quantum Monte 
Carlo-MMEH'! The HNC-EL approach has been used 
extensively in studies of the third and higher layers of 
helium films on graphite. This approach determines sta- 
ble coverages of the helium layers. Because the theory 
requires a uniform film, calculations cannot be made for 
all coverages. Breakdowns in the theory are interpreted 
as occuring at coverages where the film is unstable to the 
formation of droplet patches or to layer promotion. The 
theory predicts at least three liquid layers will form on 
top of the solid first and second layers, and yields a max- 
imum coverage value before promotion of 0.065 atom/A^ 
for each layer, in good agreement with_but somewhat 
below, the experimental value of 0.076BE£l 

In this paper we present results for the third and 
fourth helium layers using the path-integral Monte Carlo 
(PIMC) method. Our simulation is able to take into ac- 
count the effect of the second layer's corrugations and 
zero-point motion on the third layer. We are also able 
to allow for the possibility of promotion and demotion 
of particles between the second, third, and fourth layers, 
and these effects are observed. Finally, our simulation 
method can be applied to the entire range of possible 
phases in a layer, from liquid droplets to full solidifica- 
tion. Thus we are able to probe both the low density and 
high density phases of a layer. 



where p is the density matrix, is a configuration of 
N particles, and r = 13/M . The sum over P is the sum 
of all possible permutations of particle labels. Both the 
configurations and the permutations are sampled by our 
PIMC method. The number of intermediate configura- 
tions used, M , is referred to as the number of inverse- 
temperature slices. 

To implement this method, two important ingredients 
are needed. First, a starting approximation for the den- 
sity matrix at r is required. The simplest of these is the 
so-called semiclassical approximation. The drawback to 
this approach is that very high starting temperatures are 
needed to obtain an accurate approximation. Methods 
for improving the starting approxiniation can lower the 
starting temperature to — 40if .cZl 

The second important feature of the method is mul- 
tilevel sampling. The path- integral, Eq. can be 
thought of as a system of N ring polymers, each with 
M beads. Particle permutations correspond to splicing 
the ring polymers together. In multilevel sampling, we 
update particle positions over a section of the polymer 
chain. This allows us to implement permutations by 
spreading them over several inverse-temperature slices. 
The number of slices to be updated is 2', where I is the 
level of the move. The value of I must be chosen to bal- 
ance the acceptance of particle moves and particle per- 
mutations. Increasing the value of I increases the rate of 
accepting permutations but decreases the acceptance of 
new particle positions, while decreasing I has the opposite 
effect. Typically, we take I = 3 since this gives the op- 
timum balance between the acceptance of permutations 
and new particle positions. 



II. SIMULATION METHOD 



B. Third and higher layers 



Our simulation is performed using a path-integral 
Monte Carlo method that includes particle permutations 
and substrate effects. The general method for bulk sim- 
ulation has been reviewed elsewhereJH3 and our modi- 
fications for simulatijig layered systems are given in a 
previous publication.a Some of the present calculations 
have required minor changes to this method, so we will 
briefly outline the procedure now in order to explain the 
modifications. 



A. Overview of procedure 

The partition function Z for a system of N bosons at 
the inverse temperature /3 can be expanded as a path- 
integral by inserting M intermediate configurations: 
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xp(Ri, R2; t)p(R2, R3; r) . . . p(Rm, PRi; r), (1) 



In our calculations for the third and higher layers, we 
use a simulation cell that is designed to accommodate the 
first layer helium solid. Periodic boundary conditions are 
applied in the plane of the substrate. We approximate 
the effects of the solid first layer on the second and higher 
layer atoms by placing frozen atoms at triangular lattice 
sites on the substrate. These frozen atoms are located 
2.%A above the graphite surface, the heisfait of the first 
layer as indicated by neutron scattering.Es The second 
layer also solidifies before proniotion to the third layer 
occurs. In a previous simulation^ we determined that the 
highest second layer coverage before layer promotion was 
0.2117 atom/A^. This simulation was performed with 
20 second layer atoms above 30 frozen first layer atoms 
in a cell of dimensions 15.075A x 15.67A. The majority 
of our simulation results for the third and higher layers 
were obtained using a simulation cell of this size. 

In order to investigate the effects of an inert versus 
active second layer, we have performed two sets of cal- 
culations. In the first, we equilibrate the second layer at 
0.2117 atom/A^. The positions of these atoms are then 
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frozen and are no longer included in the sampling. Ad- 
ditional atoms are then placed above this inert substrate 
and have their positions and permutations sampled. We 
refer to these additional atoms as being "active". The 
bisection level used in this simulation was ^ = 3. All er- 
rors that we report for these calculations are statistical 
errors arising from the Monte Carlo simulation of the ac- 
tive atoms above a particular frozen second layer. There 
will be additional systematic errors that arise from our 
particular choice for the frozen second layer configura- 
tion. 

In the second set of calculations, both the second and 
third layer particles are included in the sampling. There 
is a potential problem with doing Monte Carlo calcu- 
lations on such a system. The third layer is liquid and 
permutations will occur at low temperatures. This favors 
using ^ = 3 for this layer. The second layer, on the other 
hand, is solid and furthermore increases in coverage as 
the overall coverage is increased. Using I = 3 produces a 
low acceptance rate for particle moves in the compressed 
second layer. Sampling efficiency can be improved by us- 
ing Z = 2 for this layer. Thus a single value of I for the 
entire system is not optimum. 

To give each layer the best value of ^, we have parti- 
tioned atoms into "second layer" and "third layer" atoms. 
"Second layer" atoms have their positions initially taken 
from an equilibrated second layer solid. "Third layer" 
atoms are started from an initial configuration of atoms 
placed at various heights above the second layer. Promo- 
tion and demotion between the layers are allowed, but the 
layer label of the atoms does not change. "Third layer" 
atoms are sampled with l^rd — 3, while "second layer" 
atoms are sampled at l2nd — 2. Permutations are allowed 
between atoms with the same layer label, but we do not 
allow atoms with different layer labels to permute. This 
limitation is not a problem for atoms demoted to the sec- 
ond layer after being initially placed on the third layer, 
since exchanges are uncommon in the solid second layer. 
The promotion of a particle from the second layer to the 
third occasionally occurs after a particularly long simu- 
lation. At most, only one particle is promoted, so this 
will have little effect on the permutations in the third 
layer in all but the lowest coverages. This approach gives 
a higher acceptance rate for second layer particle moves, 
while allowing permutations to occur reasonably often in 
the third layer. We have found that this partitioning of 
I lowers the energy by an amount ranging from 0.0 K 
to 0.1 K per atom relative to energy calculations with a 
single value I = 3 for the entire system. The smallest en- 
ergy shift occurred at the lowest coverage tested, 0.2286 
atom/A^. At this coverage the two values were within 
error bars. The greatest energy shift occurred near the 
equilibrium coverage of the third layer liquid. This was 
the highest coverage we tested. We also compared the 
energy values at selected coverages to calculations per- 
formed with a single value I = 2 for the entire system 
and found them to agree, but third layer superfluidity 
was suppressed, as expected. 



III. RESULTS FOR THE THIRD AND FOURTH 
LAYERS 

Before presenting our simulation results, we wish to 
clarify the convention that we use to report adsorbed 
helium coverages. Normally, the values we give for the 
density are for the total adsorbed helium (first two solid 
layers plus any additional coverages for the higher layers). 
Relative coverages within a layer are prefaced with a ref- 
erence to that layer. For example, the coverage 0.2583 
atom/A^ using our standard simulation cell corresponds 
to a third layer coverage of 0.0466 atom/A^ plus two solid 
layers with a combined coverage of 0.2117 atom/A^. 

A. Results with the inert second layer 

1. Layer promotion and demotion 

In order to study multiple layers of the helium film, it 
is first necessary to establish that our simulation does, 
in fact, produce distinct layers with increasing cover- 
age. This is illustrated in Fig. |^, which shows the growth 
of the density profile perpendicular to the substrate for 
the third and higher layers with increasing density. The 
peaks associated with the third and fourth layers can be 
clearly observed, as can the beginning of the fifth layer 
peak at the highest simulated coverage. Also as can be 
seen in the figure, the continued growth of the third layer 
peak for all coverages indicates compression of this layer 
even after atoms are promoted to the higher layers. 

Promotion to the fourth layer occurs for coverages 
greater than 0.2837 atom/A^ (17 active atoms). This 
coverage may be determined from the following obser- 
vation, which is illustrated in Fig. For coverages at 
and just below promotion, increasing the coverage in- 
creases the peak height but does not appreciably change 
the profile's width. Just above this density, the peak 
height docs not change, but an abrupt increase in the 
width is observed. With increasing coverage, the tail of 
the profile just above layer promotion evolves into the 
fourth layer peak. This is also illustrated in the figure. 
This value for layer promotion is in agreeiafint with ex- 
periment. Heat capacity measurementstilu show layer 
promotion at 0.288 atem/A^, while isothermal compress- 
ibility measurementgl3 give a somewhat lower value 0.280 
atom/A^ for the promotion density. 

Layer promotion is also signaled by a change in the 
density dependence of the chemical potential. By differ- 
encing our calculated total energy values, we can obtain 
the chemical potential ii. This is plotted in Fig. ^. The 
values for the energy per particle are given in Table |. 
As can be seen from the figure, the chemical potential 
increases rapidly just below layer promotion until it be- 
comes favorable to promote an atom to the next layer. 
Above layer promotion, the chemical potential remains 
roughly constant with increasing density, as it should for 
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FIG. 1. Density profiles at 400 mK for the second, third, 
and fourth layers as a function of height above the graphite 
substrate. The leftmost peak is for the frozen, equilibrated 
second layer. The coverages shown begin at 0.2329 atom/A^ 
and increase in increments of 0.085, up to 0.3345. The density 
profile for 0.3811 atom/A'^ is also shown. The profiles are 
normalized so that integration gives the number of atoms. 
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FIG. 2. Density profiles near layer promotion. The cover- 
ages are 0.2794, 0.2837, 0.2879, 0.2964, and 0.3005 atom/1^ 



TABLE I. Energy /particle versus coverage at 400 mK for 
the third layer when the second layer is frozen. All calcula- 
tions use the 15.075A x 15.667^4 simulation cell. The number 
in parenthesis gives the error in the last two digits. 



the liquid-gas coexistence region. This change in the den- 
sity dependence of /i, together with the density profiles, 
determines the promotion density. 

Layer promotion in general is related to the chemical 
potential. As the density of a particular layer increases 
above the equilibrium coverage, so does fi. Each atom 
added to the system must choose whether to go in the 
dense layer or sit above it in the unfilled layer. This 
choice is governed by the chemical potentials for the two 
layers. Above the equilibrium coverage of a layer but be- 
low promotion, the /i of the dense layer is lower than /x 
for the unoccupied layer. However, the chemical poten- 
tial for the dense layer increases rapidly with increasing 
density and at some coverage, the chemical potential for 
adding an atom to the dense layer exceeds the chemi- 
cal potential for adding it to the unoccupied layer. The 
above discussion is illustrated for the third and fourth 
layers by Fig. ^. As can be seen in the range 0.26 to 0.28 
atom/A^ each additional atom causes a rapid increase in 
II. Around 0.284 atom/A^, fi^ of the third layer exceeds 
/i4 of the fourth layer, and so the next atom added to 
the system will be preferentially promoted to the fourth 
layer. 

The beginning of particle promotion to a new layer 
does not signal the end of the filling of the old layer. As 
can be seen in Fig. ||, adding atoms to the system at 
the densities 0.306 and 0.334 atom/A^, marked by the 
letter "D" and the dark arrowheads, causes the chemical 
potential to increase above the chemical potential at layer 
promotion. Immediately above both of these densities, 
fi drops back to values equal to fi at promotion. This 
signals particle demotion to the third layer. Figure ^ 
illustrates this for densities around 0.306. Notably, each 
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additional demoted atom increases the energy required 
to further increase the layer density, so /x at promotion is 
only a rough estimate of the chemical potential required 
to demote additional atoms. 

Particle demotion may be understood as a balance be- 
tween many factors. Atoms are initially promoted above 
a layer because this is energetically favorable. The first 
promoted atom loses the large energy benefit for being 
close to the substrate, but gains kinetic energy since it 
is free to move about on the open surface, with its wave 
function no longer constrained by the other atoms. It also 
retains some of the potential energy advantage gained 
from having neighboring helium atoms. Adding more 
atoms to the system will favor the formation of droplets. 
However, once droplets have formed in the new layer, the 
next atom added to the system faces a different choice 
than the first promoted atom. This is illustrated in Fig. 
^. If it goes into the less dense outer layer, it gains some 
attraction from the other atoms in this layer. However, 
it no longer gains the kinetic energy advantage that the 
first promoted atom had. On the other hand, if the added 
atom goes into the dense lower layer, it regains the benefit 
of being closer to the strongly attractive substrate. Fur- 
thermore, it has the energy advantage for having more 
helium neighbors: those in its layer, those in the layer 
above, and those in the layer below (not shown). In con- 
trast, if the added atom goes into the outer layer, it has 
less neighbors. It should be noted that each demotion to 
the dense layer does significantly increase the chemical 
potential for this layer, so at some point it will be more 
favorable for atoms to be added to the outer layer again. 




FIG. 4. Density profiles for the third and fourth layers. 
Coverages are given in A~^. 



S 
o 

Q. 




(A) 



0.26 0.28 
Coverage(A"^) 

FIG. 3. Chemical potential for third and fourth layers. 
The arrows indicate the third layer promotion density (P) 
and the fourth layer equilibrium coverage (E). The densities 
denoted by "D" and the dark arrows are discussed in the text. 
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FIG. 5. An illustration of the circumstances of particle 
demotion. The shaded atom is added to the layered system. 
In (A), it joins the less dense outer layer. In (B), it is demoted 
to the more dense underlayer. 
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2. Third and fourth layer phases 

We expect four principal regions in the third layer be- 
fore fourth layer promotion. These are a low density gas 
phase (which will have a negligible density at low temper- 
atures), a droplet region, a metastable liquid region, and 
an equilibrium liquid phase. The droplet region consists 
of a liquid phase separated from a low-density gas by an 
interface. In the metastable region, the droplet phase is 
replaced by a stretched uniform phase which has nega- 
tive spreading pressure. The crossover from the droplet 
region to the stretched liquid phase occurs at the spin- 
odal point. Direct evidence that the layer is liquid comes 
from torsional oscillator measurements, which detect su- 
perfluidity up to layer completion. The isothermal com- 
pressiMity has been measured for the third and higher 
layersEj and exhibits a divergence that is associated with 
the spinodal point. The equilibrium liquid covpsage can 
be inferred from heat capacity measurements. til Below, 
we present evidence of each of these phases using several 
different observables. 

First, we can establish the existence of droplets and a 
uniform liquid phase at different densities with the ra- 
dial distribution function, g(r), which provides a direct 
probe of short and long range behavior. Calculations for 
the third layer are shown in Fig. ^. These are plotted 
as functions of the magnitude of the distance vector be- 
tween pairs of atoms, projected onto the plane of the sub- 
strate. These calculations of the averaged g{r) smooth 
out possible anisotropics induced by the corrugations of 
the underlying solid helium layer. The g(r) for the three 
coverages shown in Fig. ^ are representative of the droplet 
region, the equilibrium liquid, and the liquid near layer 
promotion. At the lowest coverage, 0.233 atom/A^ (5 
active atoms), the radial distribution function drops be- 
low unity at large distances, as would be expected for a 
droplet phase. The actual dimensions of the droplet for 
a given density depend on the size of the simulation cell. 
For the intermediate coverage, 0.2624 atom/A^ (12 active 
atoms), the first peak has changed only slightly, but the 
long range behavior is noticeably different, rising again 
past unity instead of dropping continuously. At the high- 
est coverage, 0.284 atom/A^, the system shows evidence 
of increased correlation, but the long range behavior can- 
not be determined due to the small size of the simulation 
cell. This largest coverage also has a different short range 
behavior, showing an increased probability that the pro- 
jected distance between two atoms will be less than 2.0 
A. In part, this is a result of the thickening of the layer, 
as can be seen in Fig. |l|. Because the atoms can be at 
different heights above the substrate, the projected dis- 
tance between them can become smaller than would be 
possible in strictly two-dimensional calculations. 

We can also gain some insight into the layer phases by 
examining contour plots of the probability distribution 
of atoms in the plane of the substrate. Plots near equi- 
librium and layer promotion are shown in Fig. 0. The 
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FIG. 6. The radial distribution function for the third layer 
at the indicated coverages, in atom/A'^. 



high density liquid shows noticeably more localization 
than the equilibrium fluid. An increased correlation can 
also be seen in the radial distribution function at high 
density. This suggests that the system may be nearing a 
liquid-solid coexistence phase. 




FIG. 7. Probability densities for the third layer liquid 
near equilibrium (0.2625 atom/A^, left) and just before layer 
promotion (0.2837 atom/A^, right). 

Having established that the third layer has gas-liquid 
and uniform liquid phases, we can next determine the 
equilibrium coverage ranges of these phases. This can be 
done by using the Maxwell construction. At low tem- 
peratures, the total energy and the total free energy are 
nearly equal, so coexistence regions may be identified by 
applying the Maxwell construction to the energy. The 
results for the low temperature (400 mK) temperature 
scans are shown in Fig. ||. We have verified that the en- 
ergy values shown are effectively zero temperature results 
by recalculating some values at 500 mK. In all cases the 
calculations at the two temperatures agreed within error 
bars. The values shown in the figure have been shifted 
by the amount NactSo, where cq = — 15.897±0.024 is the 
minimum energy per particle, and Nact is the number 
of active atoms in the simulation. At low temperatures. 



6 



the gas phase will have zero coverage and thus zero to- 
tal energy. We can thus draw a coexistence line between 
the beginning of the third layer, 0.2117 atom/A^, and 
the coverage with minimum energy per particle. This 
higher coverage is the equilibrium liquid density. We 
find this coverage to be 0.2625 atom/A^ {Nact — 12). 
The best chi-squared parabolic fit around this minimum 
gives 0.2645(9) atom/A^. The number in parenthesis is 
the error in the last digit. At this density the layer is com- 
pletely covered by a uniform liquid. Below this value, the 
system enters the gas-liquid coexistence region. The en- 
ergy values in the coexistence region lie above the coexis- 
tence line, either because the liquid phase is unphysically 
uniform, or else because of the appreciable cost for creat- 
ing a phase boundary in a finite-sized system. The third 
layer equilibrium coverage, 0.0528(9) atom/A^, is compa- 
rable to (and slightly higher than) the equilibrium|-C|Qver- 
age found for the second layer, 0.0480(6) atom/A^.&cl For 
both layers, simulated with the same size cell, the energy 
minimum occurs when the system contains 12 atoms. 

The equilibrium density that we determine is in good 
agreement with both heat capacity and torsional osciila-. 
tor measurements. In the measurements of Grey wall Jlil'El 
the low temperature heat capacity depends linearly on 
density from the beginning of the third layer to 0.260 
atom/A^, 'piiis linear dependence is a signal of phase, 
coexistence. c3 The torsional oscillator measurementsH 
provide evidence of a similar region. The temperature 
of the oscillator's dissipation peak, which gives a rough 
estimate for the superfluid transition temperature, is in- 
dependent of coverage from 0.22 to 0.26 atom/A^. This, 
is characteristic of a surface covered by liquid droplets.til 
Increasing the coverage in the droplet region increases 
the size of the droplets, not their density, and so the 
transition temperature remains constant. 

For densities above the equilibrium coverage and below 
layer promotion, we have not been able to rule out the on- 
set of solid-liquid coexistence in the third layer. Our pro- 
cedure for identifying such regions in the second layer was 
to look for the instabilities in the total energy that signal 
phase coexistence. This was possible because that layer 
solidified before layer promotion. In the present case, 
though, the energy cannot be used to find solid-liquid 
coexistence because any instabilities that might signal 
solidification are inextricably entangled with the liquid- 
vapor phase coexistence that occurs in the fourth layer at 
the same coverages. Specifically, layer promotion begins 
at the third layer coverage of 0.0720 atom/A^. The next 
increment in cove rag e that we can simulate using the cell 
described in Sec. is 0.0762 atom/A^. In our prevums 
simulation of the second layer using the same sized cellj3'Ll 
we determined that an incommensurate solid begins to 
form at the second layer coverage 0.0762 atom/A^. Ap- 
parently, promotion preempts solidification. Also of note 
is that the density profile of the third layer at 0.0720 is 
much less peaked than the second layer at the same layer 
coverage, so any third layer solid phase must have very 
large zero-point motion. One possible consequence of 



the third layer entering solid-liquid coexistence is that it 
may fully solidify under compression of further adsorbed 
layers. As we have discussed in the previous section, the 
third layer density continues to increase even after fourth 
layer promotion. Layer promotion is not a phase tran- 
sition and can occur whenever the chemical potential of 
the system favors it. Thus it is possible to have layer 
promotion in the middle of a phase transition. The third 
layer solid-liquid phase transition can be completed after 
fourth layer promotion as additional atoms added to the 
system get demoted to the third layer. We note finally 
that there is some experimental evidence suggesting that 
third layer solidification may occur. In the heat capacity 
results presented by Greywall,B a small peak at about 
1.8 K can be observed for coverages beginning at 0.3100 
atom/A^, between the rounded heat capacity feature as- 
sociated with the fourth layer liquid and the sharp peak 
associated with the melting of the second layer solid. 

The fourth layer also exhibits a self-bound liquid cov- 
erage. We can identify this in the same manner as be- 
fore. Since promotion to the fourth layer occurs above 
0.2837 atom/A^, we can consider this coverage to cor- 
respond to a zero density fourth layer gas. We apply 
the Maxwell construction again to the fourth layer and 
determine a bound hquid phase at 0.3345 atom/A^ (29 
active atoms), giving a fourth layer equilibrium density 
of 0.0508 atom/A2. The soHd line in the fi gure gives the 
maximum possible range for gas-liquid coexistence in the 
fourth layer. All energy values in this region are on or 
above the coexistence line. This coexistence region agrees 
with heat capacity measurements, cl which exhibit linear 
isotherms in the fourth layer up to 0.3300 atom/A^ at 
low temperatures. 

A trend can be noticed in the equilibrium density as 
one progresses from the inner to outer layers. For the first 
layer of helium adsorbed on a flat substrate, we have cal- 
culated the equilibrium density to be 0.0450(6) atom/A^, 
close to the two-dimensional value.t3 Successive increases 
are observed in the second and third layers, as the zero- 
point motion of the layers perpendicular to the plane of 
the substrate becomes larger. From a two-dimensional 
point of view, this motion has the effect of softening the 
hard cores of the helium atoms. 

In Fig. ^, we do not plot energy values for the cov- 
erages between 0.2413 and 0.2583 atom/A^ at 400 mK. 
Calculations that we performed for these densities at 500 
mK actually showed an energy decrease, outside of error 
bars. We attribute this to variations inherent in using 
frozen configurations for the second layer: the two calcu- 
lations required two different configurations, and so there 
will be some systematic difference in the energy values. 

We have repeated the third layer calculations at 400 
mK using a larger simulation cell (21.105A x 20.889A). 
This will allow us to examine finite-size effects and the 
effects of different configurations. We will also be able 
to examine the intermediate region that we excluded in 
the calculations using the smaller cell. The total energy 
values are given in Fig. i These energy values have been 
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N 




E/N (K) 


7 


0.2291 


-15.199(51) 


8 


0.2313 


-15.362(47) 


12 


0.2404 


-15.442(33) 


15 


0.2472 


-15.689(36) 


18 


0.2540 


-15.885(36) 


24 


0.2676 


-16.005(27) 


26 


0.2722 


-15.908(25) 


28 


0.2767 


-15.833(30) 



TABLE II. Energy/particle versus coverage at 400 mK 
above the frozen second layer using the 21.015A x 20.889A 
simulation cell described in the text. The number in paren- 
thesis gives the error in the last two digits. 
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FIG. 9. Total energy for the third layer using the larger 
simulation cell described in the text. The solid curve is a fit to 
Eq. (0) . The arrow indicates the equilibrium liquid coverage. 



70.0 



50.0 




30.0 



10.0 - 



0.22 



0.24 



0.26 0.28 0.30 
density (A'^) 



0.32 



0.34 



FIG. 8. Shifted total energy for the third and fourth layers. 
The horizontal dashed line is the gas-liquid coexistence line 
for the third layer. The solid line is the gas-liquid coexistence 
region for the fourth layer. The vertical dashed line gives the 
density of layer promotion. Gas-liquid (L-l-G) and uniform 
liquid (L) coverage ranges are indicated for both layers. 



shifted by subtracting out the gas-liquid coexistence line, 
as described above. The energy per particl e value s used 
to construct this figure are given in Table III A 2. The 
minimum energy per particle is — 16.005±0.027 K, which 
occurs at 0.2676 atom/A^ (24 atoms). We fit these en- 
ergy values to a polynomial in the form 



E/V = eo + B{ 



Po 



p- 



Po 



(2) 



where p is the coverage. Fitted parameters are given in 
Table [II. Not all digits are significant. Notable from the 



table is that the equilibrium liquid coverage occurs at the 
same density as in previous calculation with the smaller 
cell, despite the fact that there is a small discrepancy in 
the energy values at similar coverages between the cells. 
In our previous calculations on the second layer, we ob- 
served that finite-size effects on the energy were negligi- 
ble. Thus we believe that the difference in the minimum 
energy values obtained with the two cells is attributable 
to the particular frozen second layer configurations that 
we used. We also observed that a calculation using the 
larger cell but at 500 mK and 0.2540 atom/A^ (18 active 
atoms) showed an energy increase over the value obtained 
at 400 mK. 

By differencing the (unshifted) total energy values of 
the larger cell, we can obtain the chemical potential /i. 
Values obtained from both the fit and the actual energy 
values are shown in Fig. |l^. The size of the chemical po- 
tential near the minimum (— 16.88i^ ± 0.28) is in general 
agreement with the value range of -^-ito -16 K obtained 
for the third layer by Clements, et a/.EZI As noted by those 
authors, their values are sensitive to the attractiveness of 
their model substrate. 

We can determine the spinodal point by taking the sec- 
ond derivative of Eq. (0) . The isothermal compressibility 
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P 3,r 3,1116 1 6r 




eo(K) 

B(K) 
C(K) 


0.0097 ±0.51 
0.2645 ± 0.0005 
2057.3 ± 237.0 
13885.7 ± 1759.0 
1.09 



TABLE III. Fitted parameters for the polynomial fit to 
the shifted total energy values. 



is given by kt — p^^ [d ii j d p)^^ , where p is the coverage. 
This diverges when the derivative of the chemical po- 
tential is zero. Below this coverage, d^/dp is negative 
and the speed of sound becomes imaginary. This is the 
spinodal point. From the fit, we determine that this oc- 
curs at 0.2518 atom/A^, a third layer coverage of 0.0386 
atom/A^. This is in agreement with the HNC-EL cal- 
culations for the third layerO Experimentally, the spin- 
odal point can be determined from the divergence of the 
isothermal compressibility. In the third layer, this occitfs, 
near the third layer coverage of 0.03 atom/12 ^ 2 K.El 
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FIG. 10. Chemical potential calculated from the total en- 
ergy values and the fit shown in Fig. ^. 

The droplet region and the uniform coverage region 
for this simulation cell can be directly identified by the 
short wavelength behavior of the static structure func- 
tion, S(k). Figure pi] shows S(k) for the coverages 0.2472 
(below the spinodal point) and 0.2676 atom/A^ (near 
equilibrium). The upward swing of S{k) for the lower 
density for small values of k indicates the presence of a 
droplet. Near the equilibrium density, the droplet has 
entirely covered the substrate and so S{k 0) ^ 0. 

At the temperature of our simulation, superfluidity can 
be observed in the third and fourth layers. Supcrfluid 
density values, ranging from intermediate third layer to 



2.0 




0.0 1.0 2.0 3.0 4.0 5.0 

k(01)(A-') 



FIG. 11. Static structure factor at 0.2472 (circles) and 
0.2676 atom/A^ (squares) at T=400 mK. 

low fourth layer coverages, are shown in Fig. ^ The val- 
ues shown are the ratios of the superfluid density to the 
third layer coverage. A slight suppression in superfluidity 
occurs for coverages just below and just above the pro- 
motion density to the fourth layer, although the noise of 
our calculations obscures this effect. Suppression of su- 
perfluidity has been observed at similar temperatures and 
coverage values in both thci-torsional oscillator measure- 
ments of Crowell and ReppylJ'El and a simulation of helium 
on hydrogen.cB Suppression just after fourth layer promo- 
tion is a consequence of gas-liquid coexistence. In the low 
coverage region, the fourth layer consists of droplets that 
lack the connectivity to exhibit superfluidity until they 
have percolated across the surfaccEflo The suppression 
of superfluidity before promotion to the fourth layer can 
be observed more clearly when we allow the second layer 
to respond to the third layer film, as described in the 
next section. 



B. Results with active second layer 

1. Layer promotion and demotion 

In the previous calculations we assumed that the solid 
second layer could be treated as inert; that is, the re- 
sponse of the second layer to the third was ignored. It 
is known experimentally, however, that the second layer 
solid continues to be compressed after third layer promo- 
tion. Evidence of this can be found in the heat capacity 
peak associated with second layer melting.B This peak 
continues to increase in temperature for the entire range 
of the third layer, indicating an inccease in the second 
layer density. Neutron scattering^I'LlI also detects com- 
pression of the second layer by the third, with an abrupt 
restructuring at intermediate third layer coverages. Sec- 
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FIG. 12. Superfluid density values at 400 mK. The frac- 
tion is relative to the number of active atoms. The arrow 
indicates the density at which fourth layer promotion begins. 



ond layer compression is not unexpected, since the com- 
pression of the iicst layer by the growth of the second is 
well established.E3 To allow for the effects of second layer 
compression, as well as the response of this layer to the 
growth of the third, we have performed calculations that 
include both second and third layer atoms in the Monte 
Carlo sampling. 

Density profiles illustrate the compression of the sec- 
ond layer in our simulation. These are shown in Fig. |l^ 
for selected coverages. The left and right peaks in the fig- 
ure are for the second and third layers, respectively. As 
can be seen, the height of the second layer peak grows 
between the coverage 0.2296 and 0.2879 atom/A^ The 
peak height then remains constant up to and including 
the coverage 0.2964, where promotion to the fourth layer 
is visible. Finally, the second layer peak increases again 
at the highest coverage examined, 0.3049 atom/A^. 

By integrating the profiles up to the minimum be- 
tween the two peaks (at approximately 7.5 A), we ob- 
tain the second layer coverages of 0.0845, 0.0885, 0.0894, 
0.0890, and 0.0929 atom/A^ for the coverages shown in 
the figure. Thus at the lowest third layer densities, no 
demotion occurs. Beginning at intermediate coverages 
and up to layer promotion, a single atom is demoted. 
At low fourth layer coverages (the highest coverage ex- 
amined) two atoms are demoted. As a consequence of 
demotion, promotion to the fourth layer is pushed to a 
higher overall coverage. As is illustrated in Fig. |l^, we do 
not observe promotion to the fourth layer until the den- 
sity has exceeded 0.2879 atom/A^, compared to 0.2837 
atom/A^ when the second layer was frozen. This value 
will be further established below by examining the den- 
sity dependence of the chemical potential|-. This value is 
still in agreement heat capacityu and isother- 

mal compressibilityE3El measurements. We note that Zi- 




FIG. 13. Density profiles for the active second and third 
layers as a function of height above the substrate. The cov- 
erage values are 0.2296, 0.2625, 0.2879, 0.2964, and 0.3045 
atom/A^. 

The layer promotion illustrated above is accompanied 
by a discontinuity in the chemical potential. We can 
obtain this quantity by differ encing our calculated total 
energy values, given in Table III B 1 , and plot the results 
in Fig. |lj. Included in this figure are results pbtained for 
the second layer from a previous calculationB As can be 
seen, there is a distinctive change in the density depen- 
dence of the chemical potential /x around the promotion 
density. Below the promotion density, the energy changes 
very rapidly with increasing coverage. Near promotion, 
an added atom has the choice of going to the unoccupied 
third layer or the dense second layer, and will choose the 
layer that is energetically favorable. When the promotion 
coverage is reached, the chemical potential for adding the 
atom to the second layer exceeds that for adding it to the 
third layer, and so the atom is added to the unoccupied 
layer. It can be seen from the figure that the change in 
fi associated with the promotion is quite large. In this 
case, the second layer is solid and relatively dense, so we 
expect an energy gap to be associated with promotion. 
As can be seen in the density profiles (see Fig. |l^) for 
this layer, there is a range of heights above the second 
layer that is forbidden for the third layer atom. This is in 
contrast to what happens in the case of fourth layer pro- 
motion (see Fig. |^), in which there is significantly more 
overlap between the third and fourth layers. Each addi- 
tional atom added to the system will also go the outer 
layer, but the rate of the energy change will decrease, 
since these atoms are attracted to each other and form 
a droplet. Thus we see that the chemical potential just 
after layer promotion decreases. 
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2. Layer phases 
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FIG. 14. Chemical potential for the second and third lay 
ers. Third layer and fourth layer promotion densities are in 
dicated by the arrows. 



We can determine phase boundaries for the third layer 
by again using the Maxwell construction. Figure [Tsl gives 
the total energy values for the three layer system, with 
both second and third layers active . The energy val- 
ues per particle are given in Table III B 1 . As before, 
a coexistence line can be drawn between the beginning 
of the third layer, 0.2117 atom/A^ and the equilibrium 
liquid coverage. The upper endpoint of the coexistence 
region is 0.2667 atom/A^. This is the equilibrium cov- 
erage for the third layer liquid. The intermediate en- 
ergy values are in gas-liquid coexistence. Following our 
usual procedure, we have subtracted out the gas-liquid 
coexistence line from all the energy values, using the val- 
ues -30.111 ± 0.019 K and -23.910 ± 0.013 K for the 
beginning (0.2117 atom/A^, 20 atoms) and end (0.2667 
atom/A^, 33 atoms) of the coexistence region. 

As with the results found using the frozen second layer, 
we can fit the energy values to a polynomial in the form 
of Eq. (||). This is the solid curve in Fig. |l5|. Energy 
values for the very low and high coverages of the layer 
were not included in the fit. The equilibrium coverage 
determined is 0.2653 ±0.0005 atom/A^, and the spinodal 
point is found to be near the third layer coverage of 0.040 
atom/A^. 



N 




E/N (K) 


19 


0.2075 


-30.687(19) 


20 


0.2117 


-30.111(13) 


21 


0.2159 


-29.257(39) 


24 


0.2286 


-27.391(20) 


26 


0.2371 


-26.354(12) 


27 


0.2413 


-25.937(12) 


28 


0.2456 


-25.516(13) 


31 


0.2583 


-24.462(21) 


32 


0.2625 


-24.188(15) 


33 


0.2667 


-23.910(13) 


34 


0.2710 


-23.603(15) 


35 


0.2752 


-23.314(13) 


36 


0.2794 


-23.001(14) 


37 


0.2837 


-22.699(12) 


38 


0.2879 


-22.417(19) 


40 


0.2964 


-21.816(13) 



TABLE IV. Energy/particle versus coverage at 400 mK 
with positions of both second and third layer atoms sampled. 
Calculations used the 5x3 simulation cell. The number in 
parenthesis gives the error in the last two digits. 
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FIG. 15. Total energy at 400 mK for the three layer sys- 
tem when second layer atoms are included in the sampling. 
The dashed line is the gas-liquid coexistence line for the third 
layer. The arrows indicate the beginning of promotion to the 
third layer and the equilibrium third liquid coverage. Cover- 
age ranges for the gas-liquid (G-fL) and uniform liquid (L) 
phase in the third layer are shown. The solid line is the fit to 
Eq. %. 

Suppression of superfluidity before promotion to the 
fourth layer is even more prominent in these calcula- 
tions. Figure |l^ shows the superfluid coverage (instead 
of the superfluid fraction) versus the coverage for the 
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third layer. For superfluid films adsorbed on heteroge- 
neous surfaces (such as Vycor), the superfluid cpverage 
increases continuously with increasing coverage.L3 This 
does not occur in the atomically thin fluid layer we are 
simulating. Superfluidity instead increases with increas- 
ing coverage until the equilibrium density is reached, af- 
ter which there is a plateau and then a drop. The drop 
in ps is seen experimentally as well, for coverages greater 
than 0.280. 
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FIG. 16. Superfluid coverage for the third layer at 400 mK 
with both second and third layer atoms sampled. 

The reason for this suppression before layer promo- 
tion is not clear, but three related factors can contribute. 
First, increasing the density in the liquid layer produces 
excitations that should disrupt the superfluid state. This 
energy increase can be seen in Fig. || and Fig. |l^. Second, 
increased correlation between the atoms will suppress the 
permutations required for superfluidity. We observe this 
increase in the radial distribution function. Fig. ^, and 
in the probability contours. The calculated probability 
distributions also suggests the third mechanism, namely 
that the layer is partially solidifying before layer pro- 
motion. As discussed above, we have not found clear 
evidence that this occurs. 



IV. SUMMARY 

This paper has presented path-integral Monte Carlo 
calculations for the third and fourth layers of helium on 
graphite using two different treatments of the underly- 
ing solid layers. In the first approach, we treated the 
first and second layers as inert. The second layer solid 
was initially equilibrated, and then atoms were frozen 
in some particular configuration. This simplification ne- 
glects the response of the second layer to growth of the 
third layer, as well effects such as particle demotion, but 
allows us to simulate larger systems. Using this approach. 



we have produced as many as three additional layers on 
top of the solid first and second layers. The third and 
fourth layers are found to possess self-bound liquid phases 
with layer densities of 0.0528 and 0.0508 atoms/A^, re- 
spectively. These densities are in agreement with values 
that may be inferred from heat capacity and torsional 
oscillator measurements. In the third layer, below the 
equilibrium density, we determine the spinodal point to 
be 0.0386 atoms/A^. This coverage separates the unsta- 
ble region, in which droplets form on the surface, from 
the metastable "stretched" phase just before equilibrium. 
Static structure calculations confirm that droplets occur 
at low densities. The chemical potential was also calcu- 
lated, and third layer promotion was found to coincide 
with a discontinuity in the density dependence of this 
quantity. The superfluid density has been calculated for 
coverages from the third layer equilibrium to low fourth 
layer coverages. Promotion to the fourth layer is ob- 
served at the (total) coverage 0.2837 atoms/A^. The 
continued growth of the third layer is also observed after 
fourth layer promotion. The question of partial solidifica- 
tion in this layer remains open for further investigation. 
In particular, it is possible that particle promotion pre- 
empts the formation of a third layer solid, but particle 
demotions at higher coverages may increase the density 
of this layer to the point of solidification. Experimen- 
tally, there is some suggestive evidence for third layer 
solidification through compression by higher layers. We 
find that the heat capacity feature that emerges above 
0.3100 atom/A^B might be related to the melting of this 
solid. Solidification also provides a natural explanation 
for the plateaus seen in torsional oscillator measurements 
of the superfluid density of the higher layers. 

The effects of including the second layer atoms in the 
Monte Carlo updating procedure have also been exam- 
ined. These calculations incorporate the effects of zero- 
point motion in the second layer solid and allow this layer 
to respond to the third and higher layers. We observe 
particle demotion from the third to the second layer near 
equilibrium and again near the fourth layer promotion 
density. The overall coverage at which the third layer liq- 
uid phase forms increases slightly from 0.2645 to 0.2653 
atoms/ A^ when the second layer is active. The comple- 
tion density is also slightly higher, occuring at 0.2879 
atoms/ A^. A very notable effect produced by the active 
second layer was the reduction in the third layer super- 
fluid density at high coverages. This effect has also been 
seen experimentally. 
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